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Abstract. Full jet reconstruction in relativistic heavy ion collisions provides new and unique insights to the 
physics of parton energy loss. Because of the large underlying event multiplicity in A + ^ collisions, random 
and correlated fluctuations in the background can result in the reconstruction of fake jets. These fake jets 
must be identified and rejected to obtain the purest jet sample possible. A large but reducible fake rate of 
jets reconstructed using an iterative cone algorithm on HIJING events is observed. The absolute rate of 
fake jets exceeds the binary-scaled p + p jet rate below 50 GeV and is not negligible until 100 GeV. The 
variable Sjr, the sum of the jet constituent's Et perpendicular to the jet axis, is introduced to identify 
and reject fake jets at by a factor of 100 making it negligible. This variable is shown to not strongly depend 
on jet energy profiles modified by energy loss. By studying azimuthal correlations of reconstructed di-jets, 
the fake jet rate can be evaluated in data. 



1 Introduction 

Complete jet reconstruction in relativistic heavy ion col- 
lisions is essential for a direct measurement of the en- 
ergy lost by a high-p-r colored parton traversing the col- 
ored medium produced in these collisions. Measurements 
of such energy loss have been made from experiments 
at Brookhaven National Laboratory's (BNL) Relativistic 
Heavy Ion Collider (RHIC) where single particle produc- 
tion rates in y/s^ — 200 GeV Au + Au collisions are sup- 
pressed relative to binary-collision scaled p+p rates at the 
same y/s [1]. These measurements alone, however, are in- 
sufficient to constrain the details of the energy loss mecha- 
nism. For example, the current data cannot yet be used to 
determine the relative contributions of radiative and coUi- 
sional enery loss, or even determine whether jet quenching 
is a perturbative or non-perturbative process. Because of 
the large underlying event multiplicity, dN/d?7 ~ 600 for 
the top 10% Au + Au collision cross section [^, direct 
jet reconstruction at RHIC is difficult, though attempts 
are underway [3] [4]. As a proxy for (di-)jet reconstruc- 
tion, two, high-pT particle azimuthal correlations have 
been measured. Results from two-particle correlations re- 
veal that the trigger jet is essentially unmodified while 
the recoil jet is strongly suppressed [5], further lending 
support to the idea that parton energy loss is the mech- 
anism of the suppression. However, two-particle correla- 
tions have not revealed expected modifications to jet prop- 
erties [6] such as softening of the fragmentation function 
[7], broadening of the di-jet acoplanarity [5], etc. expected 
from energy loss models that reproduce the single particle 



suppression. One possible reconciliation between the sin- 
gle particle and two-particle results is that the two-particle 
results are "energy loss biased" . That is, by requiring that 
two high-pT particles be present, events are chosen where 
the jet loses little, if any, energy. This can occur when the 
jet traverses a short path length fS] or emits no gluons due 
to fluctuations in radiation (punchthrough) . 

Direct measurements of jets reconstructed via stan- 
dard algorithms employed in p + p{p) or e~^e~ collisions 
modified to handle the underlying heavy ion event, should 
reduce such energy loss biases. If the energy loss is dom- 
inated by gluon bremsstrahlung radiated predominantly 
inside the jet cone \W], the lost energy will be recovered 
by reconstructing the full jet profile. If, however, energy 
loss is non-perturbative [TT] or radiates energy outside of 
the cone [1^, a suppressed jet rate compared to binary- 
scaled will be observed. Whatever the scenario, direct 
measurements of jets will make a significant impact in un- 
derstanding parton energy loss. 

In the very near future the Large Hadron Collider 
(LHC) at CERN will deliver colhsions of Pb+Pb at Ts;;^ 
= 5.5 TeV. At this energy and nominal delivered luminos- 
ity, 20 million jets with Et > 50 GeV will be produced 
per month of Pb + Pb running 13J. These jets should be 
visible above the large underlying event multiplicity and 
can be reconstructed by jet reconstruction algorithms ap- 
propriately modified to handle the underlying event. 

The ATLAS experiment will carry out extensive mea- 
surments of jet production in heavy ion collisions using its 
nearly hermetic electromagnetic and hadronic calorime- 
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Fig. 1. Comparison of the spectrum of jets in dN/d77=2700, 
b=2 fm, HIJING events. Tlie (blue) filled circles are the ex- 
pected spectrum from binary-scaled p + p, the diamonds are 
the raw reconstructed spectrum uncorrected for jet reconstruc- 
tion efficiency and energy resolution. The dashed line indicates 
the raw fake rate evaluated from hard-cut HIJING events. The 
(red) squares indicate the irreducible fake spectrum after re- 
jection of the fake jets. 



ters and its large acceptance silicon tracking system. Tech- 
niques have been developed for iterative cone jet recon- 
struction in heavy ion collisions in ATLAS [TT. These 
techniques require subtraction of the underlying event en- 
ergy within the calorimeter prior to running the iterative 
cone jet algorithm (see Section [3]). Fluctuations in the 
background remain after subtraction and a stable cone, 
defined by the algorithm, can be obtained from random 
or correlated background. In studying any aspect of jets 
and their modification, such fake jets must be tagged and 
removed from the data sample. It is, therefore, necessary 
to study the characteristics of fake jets so as to extract a 
sample of jets "with high purity. 

This article describes the techniques that have been 
developed to identify and reject fake jets -with the ATLAS 
detector arising from fluctuations in the underlying heavy 
ion event. The article is outlined as follo"ws. In Section [5] 
the relevant aspects of the ATLAS detector are defined. 
Section [3] presents the necessary details of jet reconstruc- 
tion in heavy ion collisions in ATLAS and the identifica- 
tion of fake jets. This is followed by Section |4] where the 
results of fake jet detection are presented. A discussion of 
potential biases associated with the method of rejection 
and evaluation of rejection with data is given in Section[5l 



2 The Atlas Detector 

The ATLAS detector [15] is a large, multi-purpose, high 
energy physics experiment designed to detect rare pro- 
cesses in p + p collisions at ^/s — 14 TeV at the LHC. 
ATLAS has also been shown to be an excellent detector 
for measurements of both high px and soft observables 
in heavy ion collisions [TH]- The jet reconstruction de- 
scribed here is purely calorimetric. The ATLAS calorime- 



ter is largely divided into a liquid argon (LAr) barrel 
electromagnetic calorimeter covering |?7| < 3.2, a LAr for- 
ward hadronic calorimeter covering 3.2 < < 5 and a 
barrel tile hadronic calorimeter covering \ri\ < 1.5. Each 
calorimeter compartment is divided into longitudinal sec- 
tions providing energy readout at different depths in the 
calorimeter. The barrel electromagnetic calorimeter is di- 
vided into 3 layers. The front "strip" layer, with approxi- 
mately 6 Xo, has a typical cell (single readout) segmenta- 
tion of Ar) X A(j) of 0.003x0.1. The middle layer, with 24 
Xo, has a cell size of 0.025x0.025 in Ar] x A(f>. The barrel 
tile calorimeter has three layers with depth-dependent cell 
size; the first layer has O.lx 0.1 segmentation. For recon- 
structing jets, towers in a fixed grid of 0.1x0.1 in Ar]X A(t) 
are built from sums of cell energies in all layers. 



3 Jet reconstruction and fake jet 
identification 

To evaluate the expected jet reconstruction performance 
in heavy ion collision with the ATLAS detector, di-jet 
PYTHIA [ITj events were embedded in HIJING 18J heavy 
ion events. The HIJING events (referred to as "unmod- 
ified") were generated without jet quenching and with- 
out a maximum px cut on the hard scattering processes. 
The unmodified events result in correlated and fluctuating 
backgrounds from (mini-)jets, c — c and b — b production, 
and longitudinal string radiation all unsuppressed from 
the lack of quenching. The HIJING events containing em- 
bedded jets were propagated through the ATLAS detector 
using the GEANT4 software package, digitized to resem- 
ble raw data, and then processed by the standard ATLAS 
offline software chain. 

Jet reconstruction was performed using a seeded iter- 
ative cone jet algorithm [12] [10] with a seed tower energy 
of 5 GeV and a cone radius R = \J AcjP' + Arf — 0.4. On 
average, in unmodified HIJING events with dN/d?7=2700, 
150 GeV of background event Et is contained within the 
cone. The method used to subtract the background is to 



1. 



find regions of high energy in the calorimeter (potential 
jet regions) 

calculate the average background Et in the calorime- 
ter for each longitudinal layer i and as a function of 77, 



Erp (j, ^7), excluding the regions found in step 1 
subtract E^^^ (hv) h'om all cells 
build the 0.1x0.1 grid of calorimeter towers from these 
modified cells. 



The subtracted towers from this procedure are provided 
as input to the seeded cone jet algorithm. 

Figure [1] shows the invariant yields of jets in b = 2 fm, 
dN/d77 = 2700, HIJING events. The sofid (blue) circles 
indicates the binary-scaled p + p iet rate obtained from 
PYTHIA. The diamonds show the reconstructed spec- 
trum, including embedded PYTHIA jets and fake jets 
which survive fake jet rejection (see Section [J]). The re- 
constructed spectrum is uncorrected for efficiency, energy 
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Fig. 2. Reconstructed jets witli Et—57 GeV from PYTHIA embedded into unmodified HIJING (ieft) and from liard-cut 
HIJING (riglit). Tiie jet from tlie liard-cut HIJING is fake since no jets of this energy sliould exist in this sample (see text). 



scale, and energy resolution. The need for these correc- 
tions is clear below 70 GeV while the reconstructed and 
expected spectra agree well above 70 GeV. For these high- 
multiplicity HIJING events, 70 GeV jets have a recon- 
struction efficiency of 70%, an energy scale offset of -2.5%, 
and an energy resolution of 25%. A complete evaluation 
of the centrality and i?T-dependence of the reconstruction 
is given in Ref. [Hj. 

Fake jets were produced by generating, simulating in 
ATLAS, and performing jet reconstruction on a second set 
of HIJING events, labeled "hard-cut" . These events differ 
from the unmodified HIJING events in that the parton 
spectrum at > 100 GeV^ is not sampled. Therefore, 
jets with Et > 10 GeV are purposefully suppressed with 
the intent being that any jet with Et 3> 10 GeV is a 
fake jet. Even with a hard cut, jets above 10 GeV can 
occur from initial and/or final state radiation. What was 
also found was that longitudinal strings can radiate into 
midrapidity resulting in 10~^ jets with Et > 50 GeV per 
event. These were removed from the fake jet sample by 
excluding jets which match a HIJING parton above an 
appropriate threshold. 

The resulting reconstructed jets from the hard-cut HI- 
JING events then define the "maximal" fake rate plot- 
ted in Figure [TJ The rate of fake jets exceeds that of the 
expected jet yield at low- i^T- The ratio of fake to recon- 
structed yields is 4.2 at 50 GeV, 0.33 at 70 GeV, and 0.03 
at 90 GeV. This shows that the contribution of fake jets 
decreases dramatically with reconstructed jet Et- How- 
ever, if it is desired to measure jets at moderate Et, a 
rejection procedure is needed. 

It is possible that the fake rate is overestimated. Part of 
the hard particle production in HIJING results from lon- 
gitudinal strings, color fields between the nuclei, which, 
because of their large energy at the LHC, can fragment 



into a jet at midrapidity. Although an attempt was made 
to remove such jets from the hard-cut events, a sample of 
hard-cut HIJING events with the longitudinal string radi- 
ation suppressed was also studied. The maximal rate was 
reduced such that the ratio of the fake to reconstructed 
yield is 1.5 at 50 GeV, 0.11 at 70 GeV, and 0.01 at 90 
GeV. The longitudinal radiation, however, did not change 
the characteristics of the fake jets that were produced, nor 
their subsequent identification and rejection. 
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Fig. 3. The reconstructed jet Et dependence of SJt for fake 
jets from the unmodified HIJING sample (circles) and for re- 
constructed embedded jets in dN/dr;=2700 (squares). The bars 
indicate the root-mean-square of the distribution in each Et 
bin. The dashed line is a fit to the HIJING fake jets. A clear 
separation between the fake jets and embedded jets is seen over 
all Et- The larger SJt for the fake jets indicates that they have 
a broader energy distribution than the embedded jets. 
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Fig. 4. Tlie jet Sr-independent cut variable asjj,, the number of a from the fake jet mean (see Eqn.[31) for two different jet 
Et bins 40-60 GeV (left) and 60-80 GeV (right). There is clear separation between the embedded jets (dashed line) compared 
to the fake jets (solid line) . 
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Fig. 5. Fake jet rejection vs jet reconstruction efficiency. Each 
point represents a different cut in usjj, from rejecting crsjj, > 
-1 (lower right) to usjj, > -4 (upper left). 

4 Fake jet rejection 

Rejection of fake jets is maximized when exploiting the 
differences between the embedded and fake jets. Figure [5] 
shows a general feature of fake jets by comparing two jets 
with the same reconstructed jet Et (=57 GeV). The jet in 
the left panel is an embedded PYTHIA jet and the jet in 
the right panel is reconstructed from the hard-cut HIJING 
events. Even though the fragmentation of the PYTHIA jet 
is not uniform in 77 and (j), a clear core is seen. The fake jet 
from the hard-cut HIJING event is a cone which stabilized 
around three towers at large angle to the reconstructed 
jet axis and no core exists. Such a structure is generally 
indicative of the jets from the hard-cut HIJING events: 
little to no core and energy at large angles from the jet 
axis. 

The variable used to discriminate against fake jets is 
Sjr , which is constructed to give larger weight to energies 
at large angle with respect to the jet direction. It is defined 



^3T = ^ Ex., 

ccllGjct 



.11 sin i^ccii 



(1) 



where Et.coU is the Et with respect to the beam and i?ceii 
is the angle (in (j) and ry) of the cell with respect to the jet. 
It has a trivial jet Et dependence. The definition can be 
rearranged to yield 



SjT = {sinR)ET,i. 



jet 



(2) 



where (sini?) = Z'ceii-Er.ceU siniJceiiZ-^Tjet is the energy 
weighted average of sini?. Therefore, broad distributions 
would have a larger SjT than narrow distributions at the 
same Et - Figure [3] shows the jet Et dependence of UJt- It 
is linear for the fake jets measured from hard-cut HIJING 
events. A rising trend is also observed in the embedded 
jets, but, because jets become more collimated with in- 
creasing Et, (sini?) decreases with increasing Et, and, 
therefore, the rising trend of SJt increases more slowly 
at high Et- This has the effect that the discrimination 
between real and fake jets becomes more effective with 
increasing jet Et- 

Because of the Et dependence of SJt, the applied cut 
is an i<^T-independcnt quantity: the number of RMS from 
the fake-jet mean, defined as 



SjT ~ {^Jt) (Et) 
<j{Et) 



(3) 



where {SjT){ET) and a {Et) are the i^y-dependent mean 
and RMS, respectively, fitted to the hard-cut HIJING jets 
from Figure [21 Figure 2] shows the distribution of crsjT for 
two different jet Et bins for embedded and fake jets. By 
definition the fake jets have a mean of zero and an RMS 
of 1 while the embedded jets have an obviously smaller 
mean in this variable. The rate of fakes compared to re- 
constructed embedded jets is also seen in these plots as 
well as its strong Et dependence. 

Because of the clear separation between the embedded 
and fake jets, a cut can be applied to reject fake jets at 
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Fig. 6. Tfie Sr-dependent ratio of {^jt} in PYQUEN jets to 
embedded jets (black circles) compared to fake HIJING jets to 
embedded jets (red squares) (from Figure |3]). No difference is 
observed between the quenched and unquenched jets and it is 
much less than the difference between the fake and embedded 
jets. 



a desired level. However, since the distributions overlap, 
there is a trade-off between fake jet rejection and efficiency 
loss. For example, the cut applied in Figure [1] is to reject 
jets with (JsjT > "2.5. For both Et bins shown in Figure 
0] this cut is not only effective in rejecting nearly all fake 
jets but also results in a sizeable loss of real jets. Figure 
[5] shows the anti-correlation between the rejection of fake 
jets and the reconstruction efficiency of embedded jets. In 
the end the statistics and the desired purity of the data 
sample will dictate what cut will be applied in real data 
(see Section 



5 Results and Discussion 

The variable Sjx has been shown to be effective at reject- 
ing fake, reconstructed cone jets from HIJING events. In 
this section the effect of energy loss on the rejection and 
evaluation of the fake jet rate with real data are discussed. 



5.1 Effect of energy loss 

In radiative energy loss, the energy profile of a given jet 
is modified jet-by-jet, and on average the energy profile 
is expected to show large modification [21]. Because the 
SJt is a shape cut, it is possible that such a cut could 
bias results from measuring the modification of the energy 
profile of a jet. 

An estimate of the effect can be obtained using Equa- 
tion [2] and assuming a Gaussian distribution of energy 
within a jet. For a Gaussian with an RMS of ~0.2, then 
(sini?) ~ (R). At a fixed jet Et the difference between 
embedded and fake {Sjr) then is directly related to the 
difference between {R). From Figure[3l at 50 GeV the fake 
jet {Sjr) is twice that of the embedded jet. Therefore, (i?) 
of the fake jet is twice that of the embedded jet. This is 



quite a large change for a jet with an RMS of 0.2 to one 
with 0.4. Such a difference between vacuum and medium 
modified shape may be expected for 50 GcV gluon jets 

An attempt at studying the energy loss effect on the 
rejection of Sjt has been done by simulating and recon- 
structing a set of PYQUEN [2? jets in the ATLAS detec- 
tor. PYQUEN is a set of routines which modifies PYTHIA 
partons by energy loss and including the radiated gluons 
prior to fragmentation. The events were generated assum- 
ing b=0 fm and they were not embedded into HIJING. 
Figure [6] shows the ratio of {Sjx) as a function of Et 
{i.e. Figure [3|) for quenched to unquenched jets (circles) 
and for hard-cut HIJING jets to embedded jets (squares). 
Little change between the quenched and unquenched jets 
is observed in contrast to the large difference between the 
fake and embedded {^Jt). 

The jet energy profile in PYQUEN jets are modified, 
but not to the same extent as those predicted in Ref . [2T] , 
leading to much less change in (R). These results were pro- 
duced after a simulation through the ATLAS calorimeter. 
It is possible that the segmentation is also coarse enough 
to smear any effect generated by PYQUEN. It should be 
clear, however, that because ZJJt is a shape variable, it 
may show different sensitivities to different jet profiles. 
The final resolution of how sensitive must wait for analy- 
sis of real data. 



5.2 Rejection in data 

Determining the absolute rate of fake jets in real data, re- 
gardless of the rejection variable, can be done by studying 
the A(p distribution of the pairs of reconstructed jets. This 
is a useful variable because 1) the position resolution of 
jets is good, 2) it has a clear signal peak at A(p — tt, and 
3) fake jets will be randomly associated with other real or 
fake jets in the events and result in a pedestal. The rejec- 
tion variable distribution should be different for different 
regions of A(j). 

This is the case for asjr- Figure [7] shows the distribu- 
tion of all jets pairs with E^ > 60 GeV correlated with 
with Ei^ > 40 GeV with different crsjT cuts applied. For 
the successively hard asjr cuts the pedestal is suppressed. 
The stars indicate the cut used in Figure [1] which removes 
a significant fraction of fake jets. With the tightest cut, 
the signal peak at Acj) ~ tt becomes slightly narrower. 
This may indicate that the hardest cut removes of the 
radiated jets in a 3-jet event. 



6 Summary 

A large fake rate of cone jets from HIJING events is ob- 
served, exceeding the expected p-|-p jet rate at low Et- 
They result from stable cones being fitted to random and 
correlated background fluctuations. The variable TIjt has 
been used to identify and reject these fake jets. It is shown 
that this variable is insensitive to the effects of energy loss 
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Fig. 7. Z\(;/) distribution between reconstructed jets witfi > 
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on jets. It is conceivable that neither the extreme fluctua- 
tions nor the characteristics of fake jets (broader, large an- 
gle energy distribution) generated in HIJING are observed 
in the data. In fact, the fake rate will strongly depend on 
the model of the background and, indeed, will require mea- 
surement in the data itself. Even so, because correlated 
and random fluctuations can produce stable cone jets, it 
is necessary for any jet analysis in heavy ion collisions at 
RHIC or the LHC to demonstrate a lack of contamination 
of such jets. 
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